Coherent Light Scattering from Semicontinuous Silver Nanoshells near the 

Percolation Threshold 
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We report on measurements of visible extinction spectra of semicontinuous silver nanoshells grown 
on colloidal silica spheres. We find that thin, fractal shells below the percolation threshold exhibit 
geometrically tunable plasmon resonances. A modified scaling theory approach is used to model the 
dielectric response of such shells, which is then utilized to obtain the extinction cross section in a 
retarded Mie scattering formalism. We show that such spherical resonators support unique plasmon 
dynamics: in the visible there is a new regime of coherently driven cluster-localized plasmons, while 
crossover to homogeneous response in the infrared predicts a delocalized shell plasmon. 



Noble-metal nano-scale shells are comprised of thin 
gold or silver films surrounding sub-micron dielectric 
cores. In recent years nanometer-scale metal parti- 
cles have been the focus of extensive studies, owing 
mainly to large enhancements of surface-induced electric 
fields at the plasma resonance of the nanoparticles 0. 
This unique optical response is well exemplified in the 
nanoshells' extinction spectra, which are governed by a 
geometrically tunable plasmon resonance. When the core 
diameter is in the sub-micrometer range, the optical re- 
sponse of the composite particle is tunable over the en- 
tire visible and near infrared spectrum Q- This consti- 
tutes a powerful tool for custom-designing Raman [3j and 
surface-plasmon- based ultra sensitive optical sensors . 

The strong optical resonance observed in spherical 
nanoshells is attributed to a delocalized shell plasmon, 
typically well-modelled by classical linear response to 
electromagnetic (EM) plane- wave scattering [5|. This 
was formulated by Aden and Kerker (AK) |6j as an exten- 
sion to classical Mie scattering theory. More refined the- 
ories utilizing ab-initio quantum mechanical modelling of 
metallic nanoshells have also verified this approach [jj. 

The existence of a spherical-shell plasmon relies on a 
symmetric and uniform shell morphology, and is typically 
observed in smooth, continuous nanoshells. It has been 
previously shown that extinction spectra of incomplete, 
highly fragmented nanoshells do not exhibit geometrical 
resonances but are dominated instead by absorption 
resonances of the metal clusters on the spheres 0. 

In this Letter we present light scattering experi- 
ments of silica spheres coated with discontinuous silver 
nanoshells. We distinguish two shell morphologies with 
differing optical signatures. For thin, two-dimensional 
(2D) fractal shells, the optical response of the metal is 
best described within the framework of percolation the- 
ory. While these shells support mainly localized plasmon 
modes and are electrically insulating, distinct geometri- 
cal resonances are clearly present, due to coherent opti- 
cal driving of cluster-localized plasmons. For thick (3D), 
conducting and porous shells classical, Maxwell Garnett 
(MG) theory is used instead to obtain the effective di- 
electric response of the shell. In both cases, although 



the shells are locally highly irregular, observed fine res- 
onances are still well modelled by the AK solution to 
Maxwell's equations for uniform core/shell systems. 

Colloidal silica spheres (0.3-1/im diameter, 2% polydis- 
perse) were coated with nanocrystalline silver shells using 
a Tollen's reagent method to form dielectric core-metal 
shell particles |lOj |. Aqueous ammoniacal silver nitrate 
was mixed with an aldehyde reducing agent such as glu- 
cose, leading to rapid precipitation of silver atoms. At 
3°C the reaction kinetics are dramatically slowed, result- 
ing in controlled growth of nanocrystalline silver, only at 
the spheres' surfaces. Typical reactions last 20-30min. 
Repeated centrifugation and sonication are utilized to 
eliminate residual silver nanocrystals from solution. De- 
spite a lack of surface functionalization, the silver adheres 
strongly to the spheres and does not separate during this 
process. More than 70 such batches were synthesized. 

Transmission electron micrographs (TEM) of more 
than 200 coated spheres per data set are used to deter- 
mine average shell thickness and polydispersity. For thin 
coatings (~20 - 30nm) ensemble polydispersities can re- 
main as low as 2%, indicating highly uniform deposition. 
High resolution TEM analyses of isolated metal clusters 
of ~ lOOnm mean size indicate the latter consist of ag- 
gregated nanocrystals 5 - lOnm in diameter. 

Further characterization of the shells' morphologies 
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FIG. 1: (a) Scanning electron micrograph of lfim diameter 
silica sphere coated with 20nm-thick silver shell. Scale bar is 
500nm. (b) Digitally processed (see text) magnified view of 
region inside the square in (a). Scale bar is 200nm. Silver ap- 
pears white and dielectric inclusions are black, (c) Histogram 
of p, generated from 80 images similar to (b). 
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utilizes high-resolution scanning electron micrographs 
(SEM). These reveal the non-contiguous nature of the 
thin shells, which we characterize in terms of their metal 
filling fraction, p, and a fractal dimension, Df. SEM 
images as in Fig^a) are first processed with a high 
pass fast Fourier transform filter to eliminate nonuniform 
background intensity variations across the sphere surface. 
Sphere curvature and shadowing effects are minimized 
by utilizing only a small portion of the exposed surfaces 
close to nadir, shown enclosed in the box in Fig^a). 
The size of this area must be small enough to minimize 
shadowing effects, while still sufficiently large to provide 
statistical integrity. With a resolution-limited pixel size 
of ~ 2nm, for coated spheres of radius R a box size of 
(i? 2 /4) ~ 200 2 pixels satisfies this condition. Linear con- 
trast enhancement and thresholding of the cropped sec- 
tion provide the binary image in Fig^b). Visual in- 
spection of such images suggests a strong resemblance to 
percolating planar metal films deposited under ultra-high 
vacuum, where both linear and nonlinear optical charac- 
teristics have been previously addressed [ll|, |l2| • 

We measure the fractal dimension of thin films such as 
in Fig^b) by implementing a box-counting method [l^j ] 
over two orders of magnitude, and find Df — 1.72 ±0.06, 
indicating diffusion- limited aggregation in 2D 0] . Since 
this is below the fractal dimension of 2D percolating clus- 
ters, Dp = 1.89 0| , we conclude these shells lie slightly 
below the percolation threshold. We have verified this 
by measuring the average filling fraction p of thin shells 
utilizing the processed SEM images. For an ensemble of 
~100 spheres we find p — 0.55, shown in FigQIc). Com- 
paring to the experimental percolation threshold, p c — 
0.68 hj, our thin shells have p < p c . 

This is supported by measurements of the bulk re- 
sistance of silver-coated spheres. A dense aggregate of 
spheres is formed by drying several concentrated drops of 
the sphere suspension between two gold electrodes sepa- 
rated by ~100/im. For shells 20 - 30nm thick we measure 
a resistance i?~20Mf2, while shells of thickness >70nm 
exhibit a resistance of i?~lf2. These are typical values, 
also obtained previously for the sheet resistance of planar 
thin metal films near the percolation threshold [llj . 

To optically characterize our core/shell composites we 
measured the extinction of freshly prepared, dilute aque- 
ous suspensions using a UV/vis spectrometer. Spectra 
of sparsely coated spheres as in FigJSfa) exhibit only a 
single peak centered near 420nm. This is a signature 
of the excitation of the dipolar plasmon eigenmodes of 
individual, non-aggregated silver nanocrystals nucleated 
on the spheres |9j] . Due to the dilute metal surface- 
coverage, these modes do not couple to the weak Mie 
resonances of the dielectric spheres, and are therefore in- 
dependent of sphere size. At the opposite extreme are 
densely coated spheres as in Figl^b), exhibiting multi- 
polar core/shell resonances In this limit the extinc- 
tion is well-modelled by applying MG theory within the 
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FIG. 2: (a) Extinction of sparsely coated 1/j.m diameter 
spheres (TEM inset) showing a dipolar plasmon resonance at 
420nm. (b) Measured extinction (dots) and fitted AK model 
using MG theory (line) for 70nm-thick shells (SEM inset). 

AK framework, as discussed in a later section. 

Nucleation and aggregation of silver crystals onto the 
spheres result in an intermediate regime of discontinuous 
nanoshells of thickness 20 - 30nm. Extinction spectra of 
the spheres in Fig^a) reveal several distinct resonance 
peaks, all shifted from the single-particle plasmon reso- 
nance, as shown in Fig^a). These are core/shell multi- 
pole modes, typically observed in high-contrast layered 
spherical particles. In addition, the inhomogeneously 
broadened enhanced red extinction tail present in fractal 
metal aggregates is strongly modulated here. These 
indicate cluster-localized plasmon modes are interacting 
via cavity Mie resonances. Furthermore, varying the core 
size affects both peak number and positions, verifying the 
geometric nature of the resonances. In this new regime of 
coherently driven localized plasmons collective response 
of a strongly disordered shell is thus achieved. 

The often-used AK formalism, which accounts for mul- 
tiple scattering in layered concentric shells Q, requires 
knowledge of the dielectric functions of core, shell and 
embedding medium materials. Using tabulated values for 
the dielectric function of bulk silver to model the extinc- 
tion coefficient exhibited serious discrepancies and failed 
to reproduce the experimental data, as seen in Fig. Ela). 

As previously indicated, thin-shell samples are insulat- 
ing, with filling fractions close to the percolation thresh- 
old. Hence, we develop a method for modelling the opti- 
cal response of such quasi-shells based on scaling theory 
(ST) of percolating planar metal films. The details of ST 
are described elsewhere [llj. and we summarize only its 
germane points. Scaling properties of fractal networks 
of conductors allow the definition of a scaling function 
|l7j |. A dielectric function e(uj;L^) is obtained in terms 
of this function and a characteristic length scale L^(lo) : 
describing the largest length scale over which homoge- 
neous averaging may occur. L^{u>) is determined by the 
smallest relevant length scale: the percolation correlation 
length, £, or the optical coherence length L(u>), given by 

6 oc (\p-p c \/p c y u (la) 
L(u) oc w-VC 2 ^) (lb) 
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threshold. In the IR, Re[e s heii(w)] crosses over to negative 
values. For common embedding dielectrics the crossover 
to metallic behavior occurs at A ~ 1.5 — 2fim. This is re- 
lated to a crossover from inhomogcncous to homogeneous 
optical response, calculated to occur near ~ 650nm in 
our samples For A ;§> A^ the shell may be treated 

as optically homogeneous. The consequences of this are 
discussed in a later section. 

The physical properties of composite materials are of- 
ten treated using an effective-medium approach. Accord- 
ing to this method, when spatial inhomogeneities mani- 
fest on scales much smaller than the relevant length scale 
in the system, here the optical wavelength, average ho- 
mogeneous treatment is justified. However, close enough 
to the percolation threshold this is known to break down 
[ill HH Eof . We have examined Bruggeman effective 
medium theory (EMT), best suited for systems near per- 
colation 0,H2]> as a plausible model for the shells' re- 
sponse. The effective dielectric function, e e ff is given by 



FIG. 3: (a) Comparison of normalized, averaged EMT (blue), 
and ST (red) extinctions with measured spectra (dots). Also 
shown is the extinction obtained using the tabulated bulk 
silver dielectric function (black), (b) Normalized extinction 
spectra of two typical thin-shell samples with same core/shell 
size. Inset: Spectra after subtracting computed background 
functions, (c) Comparison of background-subtracted, aver- 
aged extinction for EMT (blue), ST (red) and a typical data 
set (black). All samples had 1/im diameter cores and 20nm- 
thick shells. 

For a 2D system near percolation the critical exponent 
is v = 4/3, and 9 = 0.79. A dielectric function e, (e c ), 
corresponding to the thin-shell insulating (conducting) 
component [ll|. is calculated for each length scale L^(u) 
using parameters relevant to our system [1£(. 

The main difference between near-percolation planar 
films and our nanoshells stems from the closed shell 
geometry, which introduces correlations between film 
components, otherwise uncorrelated in an open, planar 
configuration. This allows coherent driving of cluster- 
localized plasmon resonances through multipole electro- 
magnetic cavity-modes, clearly resolved in the extinction 
spectra. We account for this cluster-coupling by defining 
the shell-averaged dielectric function 

eshellM = =P^c(^,L^) + (l-p)ei(oj,L^) (2) 

which is then utilized to calculate the extinction cross 
section in the AK formalism. 

The resulting effective dielectric function exhibits char- 
acteristics strikingly different from those of both bulk and 
nanocrystalline silver. Most importantly, in the visible 
and near infrared (IR) spectral range Eq. J2J yields posi- 
tive values for Re[e s heii(w)] and relatively large values for 
Im[e s heii(w)], implying a lossy dielectric response. This 
is typical of disordered 2D metals below the percolation 



P — Tj \ — + I 1 - p) — 7: x — = (3) 

ge m + (1 - g)e eS ge d + (1 ~ g)e cf[ 

where e m and td are known dielectric functions of the in- 
terpenetrating metal and the dielectric, respectively, and 
p c = g = 0.68 is the depolarization factor in 2D (llj . 
Using measured values of p we obtain an expression for 
e c ff, which is then utilized to calculate the extinction cross 
section as described above. For the relevant experimental 
parameters |l8j , e c s is a positive function over the entire 
EM spectral range, implying non-metallic response. As 
a final step for both ST and EMT approaches, the cal- 
culated extinction cross sections are averaged over the 
measured filling fraction p and sphere size distributions 
[23j |. and are plotted together for comparison in Fig|3{a). 

Comparing the models reveals discrepancies in the gen- 
eral trend of the functions, with ST better following the 
data. Most importantly, the number of resonances, their 
position and especially contrast, are poorly described by 
EMT. We note data collected from different batches may 
exhibit opposite trends, as illustrated by the two typical 
examples in Fig|3{b). This is attributed mainly to var- 
ious degrees of nanocrystal aggregation present in each 
batch, where increased aggregation results in a broad, en- 
hanced red absorption tail. To account for this we first 
construct a background continuum function by applying 
a 50nm boxcar average to the extinction spectra. Sub- 
traction of each continuum function from its respective 
spectrum reveals excellent agreement between the promi- 
nent spectral features, shown in Figl^b) (inset). 

FigEI c ) shows background-subtracted spectra, which 
may be further analyzed using bandpass spectrophotome- 
try. Due to the large number of parameters used by ST as 
well as the highly complex nature of the multipolar vec- 
tor spherical harmonics, conventional multivariate fitting 
methods proved inefficient. Nevertheless, ST clearly best 
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FIG. 4: (Color online) Calculated extinction for ST (line) and 
EMT (dashed) for p — 0.55, 20nm-thick silver shells on l(j,m 
diameter silica cores. 



reproduces the experimental data, while EMT lacks the 
contrast to account for the observed resonances. We at- 
tribute the remaining discrepancies between ST and the 
data to a non-ideal experimental system, its main limi- 
tation being irregular shell boundaries, whose additional 
losses cannot be accounted for within the AK formalism. 

As shell thickness increases through deposition of ad- 
ditional silver, we observe shifting of the resonances as in 
FigEtb) . High-resolution SEM micrographs indicate the 
increase in thickness is accompanied by filling of inter- 
particle voids, as well as shell-coarsening. Shells thicker 
than 50nm are therefore better treated as 3D aggre- 
gates of nanocrystals, with filling fractions estimated at 
0-8 < p < 0.9. Average measured polydispersities of such 
shells may still be as low as 5%. Extinction spectra re- 
veal characteristics closely similar to those of continuous 
silver nanoshells. We may therefore describe the effective 
dielectric response of porous shells using the classical MG 
model. The dielectric function e MG is given now by 
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3(1 ~ ~ £mV (fj + 2e TO ) 

1 - (1 -p){e d - e m )/(e d + 2e m ) 



(4) 



describing a solid silver shell with 1 — p volume fraction 
of spherical dielectric inclusions. Following our procedure 
for calculating the size-averaged extinction cross section 
0, we obtain the plot in Fig|2fb). We find this model 
reproduces well our experimental observations. 

An interesting discrepancy between ST and EMT 
emerges when the extinction calculations are extended to 
the IR, as in Fig0J We find ST predicts broad IR absorp- 
tion, while EMT exhibits a decaying Rayleigh tail. The 
broad peak near 3.2/zm is a mostly dipolar shell plasmon 
resonance, verified by numerical mode-decomposition of 
absorption and scattering cross-sections. This is a conse- 
quence of the above mentioned ST-crossover to metallic 
behavior, which enables the excitation of a delocalized 
shell plasmon. Note this IR resonance should be distin- 
guished from anomalous IR absorption in metal-dielectric 
composites ^|| . The latter originates from enhanced ab- 
sorption in the embedding dielectric, due to intense local 
electric fields j2^|- Here we assume both silica core and 
dielectric inclusions in the shell are lossless, therefore all 
absorption is due to intrinsic losses in the metal shell. 



In summary, we have shown that extinction proper- 
ties of semicontinuous silver nanoshells grown on dielec- 
tric cores unambiguously differ from previously reported 
continuous core/shell systems. We have demonstrated 
tuning of the plasmon resonance with shell growth - 
from a single dipolar mode in sparse coatings to multiple 
high order shell excitations in densely packed thick shells. 
We have identified an intermediate regime of thin frac- 
tal shells near the percolation threshold, where geomet- 
rically tunable plasmon resonances persist despite disor- 
der and the shells' discontinuous character. A modified 
approach for modelling the dielectric response of fractal 
shells, based on scaling theory is consistent with exper- 
imental observations. We attribute this new resonant 
regime to cluster-localized plasmons, coherently driven 
by micro-resonator modes. Future work will address the 
predicted delocalized plasmon resonance, resulting from 
a crossover to homogeneous metallic response in the IR . 
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